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c§ i ABSTRACT 



We provide a set of general tools for studying the alignments of dark mat- 

lO i ter halos and galaxies with respect to the large scale structure. The statistics 

O ■ 

of the positioning of these objects is represented by a Probability Distribution 
' Function (PDF) of their Euler angles. The mathematical identities relating 

Q\ | this PDF to the alignment of the halo axes with respect to the direction that 

characterises locally the large scale structure are given. The PDF correspond- 
ing to halos located in the shells of the cosmic voids is inferred from previous 
5_i ' results. This PDF is used to show how to recover the outcomes found for the 



alignments of the axes of these halos in simulations. We also explore the ori- 
entation of the angular momentum of the halos, both with respect to the halo 
axes and with respect to the large scale structure. We present an expression 
which describes well numerical results for the alignment of the angular mo- 
mentum of the halo with respect to the halo axes for randomly chosen halos. 
Combining that expression, with the PDF of the halo positioning with respect 
to the large scale structure we find, in the case of halos in the shells of voids, 
an alignment of the angular momentum that is opposite to that found in sim- 
ulations. To solve this issue, we propose a model that relates the orientation 
of the angular momentum with the halos axes accounting for the orientation 
of the halo axes with the large scale structure. This model is shown to recover 
accurately the observed PDF of the halo angular momentum with respect to 
the void radial direction. This model also predicts a substantial dependence 
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of the intrinsic alignment of the angular momentum on environment. In ad- 
dition, we give an expression for determining the degradation of the angular 
momentum intrinsic alignment when observational errors are accounted. This 
expression is also used to determine the departure of the observed value of 
the alignment from the initial expectation (as provided by the tidal torque 
theory) due to the rotation of the angular momentum of the halo with respect 
to the initial torque. For voids, we find that the strength of the alignment is 
reduced to half the original value. We discuss how to adapt the void results 
to other cosmic large scale structures (i.e. filaments, walls, etc). 

Key words: 

dark matter - galaxies: halos - galaxies: spiral - galaxies: structure - large- 
scale structure of universe - methods: statistical 

1 INTRODUCTION 

Our current galaxy formation paradigm states that galaxies acquire their angular momen- 
tum by a mechanism known as tidal torquing (Sciama 1955; Peebles 1969; Doroshkevich 
1970; White 1984). In this paradigm, tidal fields generated from the large-scale structure 
excerpt a torquing on the protogalactic object prior to gravitational collapse. The acquired 
angular momentum is conserved during collapse and their value will dictate the formation 
and orientation of the galactic disks. One important consequence of the tidal torquing theory 
(TTT) is an expected relation between the large-scale structure and the orientation of the 
angular momentum of the galaxies within it (see a recent review by Schafer 2009). 

From the observational point of view, Navarro et al. (2004) found a notable (with a 
significance of 92%) excess of edge-on spiral galaxies (for which the orientation of the disk 
rotation axis may be determined unambiguously) highly inclined relative to the supergalactic 
plane. Trujillo et al. (2006) perform a measurement of the rotation axes of galaxies situated 
on the shells of large cosmic voids and found that the angular momentum of these objects 
tend to be aligned perpendicular to the radial direction of the void with a significance of 
99.7% (see also Slosar & White 2009). Lee & Erdogdu (2007) combined a reconstruction of 
the density tidal shear field from the 2MASS redshift survey with the spin field derived from 
the Tully galaxy catalogue and found a correlation strength at the Qa significance level. 

From the numerical perspective, the orientation of the halos within their large-scale 
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cosmic structure has been investigated by several authors (Bailin & Steinmetz 2005; Patiri 
et al. 2006; Brunino et al. 2007; Aragon-Calvo et al. 2007; Cuesta et al. 2008). In particular, 
Patiri et al. (2006) and Brunino et al. (2007) investigated the alignment of dark matter 
halos within the shells of large cosmic voids and provided evidence that the minor axis of 
the dark matter halos points preferentially to the centres of the voids. The intermediate and 
major axis were found to be preferentially oriented perpendicularly to the void direction. In 
addition, Cuesta et al. (2008) found, with high significance, a preferential alignment between 
the angular momentum of the halo and the void direction. As for real galaxies, the angular 
momentum of the halos tend also to be aligned perpendicular to the radial direction of the 
void. 

Although the orientation of the galaxies and halos within the large-scale structure has 
been identified both observationally and numerically, there has not been any previous quan- 
titative study neither of the relationship between the various kind of alignments already 
determined, nor of the connection between these alignments and the existing results for 
the positioning of the angular momentum with respect to the halo axes. In this paper we 
use and discuss several previous results concerning the alignments of halo axes and angular 
momentum with respect to their corresponding voids for those halos lying in the shells of 
the voids. Not all these results are independent, so we intend to establish a convenient set 
of independent results from which all other results may be obtained. To this end we shall 
present the needed mathematical identities relating the various kind of alignments that we 
shall consider, and use them to infer the set of independent results that we shall choose. 
Then we shall show that, using the appropriate mathematical relationships, all previous re- 
sults can be derived, thereby proving their consistency. In this way we shall not only present 
previous results in a compact and ordered fashion, but shall also present results whose form 
is valid for the alignment with respect to features of the large scale structure other than 
voids. 

In Section 2, we show how to characterise uniquely the statistics of the orientation of 
the halos in the shells of the voids with respect to their corresponding voids and how the 
alignments of the various halo axes may be derived from it. In Section 3, we provide a moti- 
vated model for the probability distribution for the orientation of the angular momentum of 
a randomly chosen halo with respect to its axes. Then, through a modification of this model, 
we obtain a model for the alignment of the angular momentum with respect to the halo axis 
for a halo in the neighbourhood of a void (or, in general, within a large scale tidal field). 
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Using this model along with the statistics for the orientation of the halo with respect to 
the void we obtain the probability distribution for the alignment of the angular momentum 
with the radial direction of the void. Then, comparing this distribution with that found in 
simulations, we calibrate that model, inferring a very substantial dependence of the align- 
ment of the angular momentum with respect to the halo axes on the environment. Finally, 
in Section 4, we discuss the relationship between the alignment of the angular momentum 
of the halo at the moment of turn around and at present, and how any uncertainty in the 
determination of the angle of the rotation axis of galaxies with respect to the radial vector 
joining the centre of the corresponding void with the galaxy diminishes the strength of the 
alignment. A discussion and a summary of conclusions are given in Sec. 5. 

2 ALIGNMENTS OF HALO AXES WITH RESPECT TO THE VOID 
RADIAL DIRECTION 

We shall first consider the alignments of the three halo axes with respect to the void radial 
direction. These alignments were quantified first by Patiri et al. (2006) and, with better 
statistics, by Brunino et al. (2007) and Cuesta et al. (2008). In this last work it was argued 
that all positional information of halos with respect to the void may be encoded by the three 
Euler angles (6, ip, 0) that transform the three axes of the halo into the local orthogonal 
base (i.e. the void radial direction, and the directions along parallels and meridians for any 
arbitrarily chosen set of geographical coordinates on the shell of the voids). 

We call 6 the angle between the major axis and the radial vector joining the centre of the 
void with the centre of the dark matter halo, ip the angle between the direction perpendicular 
to the major axis within the plane of the void and the minor axis, and <fi the angle between 
the direction perpendicular to the major axis within the plane of the void and the direction 
along the local parallel (see curve A in Fig. 1). Given the isotropy of the problem within 
the plane of the void, must be uniformly distributed between and ir. This was checked 
in simulations, for consistency, by Cuesta et al. (2008). Note that since these Euler angles 
determine the position of three orthogonal axes rather than a set of three orthogonal vectors, 
we have 0e[O,7r/2], ipe[0,7r] and 0G[O,7r]. 

If the shape of the halos were not aligned with any particular direction, cos#, ip, and 
4> will be distributed uniformly. However, the numerical simulations have shown that halos 
are not oriented randomly within the large scale structure. The probability distribution of 
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Dark Matter Halo 




Void 



Figure 1. Schematic illustration of the orientation of the dark matter halo in the void surface. The position of three axes of 
the dark matter halo are characterised by the Euler angles (8, ip, </>). A, M are the local parallel and meridian, respectively. 



9 has been explored by mean of several simulations. In the particular case of halos within 
a shell with a thickness of 4 h _1 Mpc around voids larger than 10 h _1 Mpc (defined by not 
containing halos with masses larger than 7.6 x 10 11 1i _1 Mq), the simulations have shown 
that this probability distribution can be well described by the following expression: 
pdfi 



jj, = cos9 



(1) 



[1 + {p 2 - l)/X 2 ] 3 /2 

with p a free parameter that is found to be 1.111±0.004 (see e.g. Cuesta et al. 2008). 
Cuesta et al. (2008) also have studied in the simulations the distribution of the angle if). 
They found that the distribution have a small, but highly significant, departure from unifor- 
mity. Given the small amplitude of this non-uniformity, we propose the following expression 
containing the simplest ip dependent term which is even with respect to ip — 7r /2. This is a 
condition that must be satisfied due to the symmetry of the problem. 



P(ip)dip = -(1 +/?cos(2^))dV 



7T 



(2) 



This expression is able to describe fairly well the data from the simulations with f3=- 
0.093±0.005. This means that it is more likely to find the minor axis perpendicular to the 
intersection of the plane of the void with the plane defined by the middle and minor axes 
than parallel to this intersection. The opposite is obviously true for the middle axis. 

Note that in the above expression ip is assumed to take values between and n/2. As 
we mentioned above, tj) actually goes between and ir, but due to the symmetry of P(V>) 
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around ip=n/2 it is not necessary to consider values of i[) larger than 7r/2 (even though the 
position ip is not, in general, equivalent to ir-ip). However, we would like to stress that in 
the process of deriving the forthcoming mathematical identities tp is assumed to go from 
to 7r and it is only in the final expressions that the symmetry of P(ip) is exploited, so that 
in all the final expressions ip is assumed to be between and ir/2. 

From the two expressions we have stated before, it follows that the full statistical infor- 
mation about the orientation of the halos with respect to the void is encoded in the joint 
probability distribution for cos(6 l ) and ip: P or (cos #,?/>). We drop the dependence on <fr since 
it is trivial. We propose the following expression to combine the above two probabilities, 
P(cosfl) and P(ip): 



The first parenthesis corresponds to the probability distribution for cos# (Eq. [Q and 
the second parenthesis to the conditional probability distribution P(ip\9) for ip for a given 
value of 9. To relate the parameter f3' with f3 we applied the fact that, for the halos under 
consideration (p=l.ll), marginalising the above probability P or with respect to cos#, Eq. [2] 
must be obtained. On doing this we obtain: /3'=1.483/3. 

We may now note that the same physical cause (the shear of the velocity field) originates 
the departure of the distribution of cos 9 and i[) from uniformity. Both effects are of the same 
order on the shear. Consequently, to leading order the quantity determining the departure 
from uniformity of cos 9 (i.e. p — 1) must be proportional to the quantity determining the 
departure of if) from uniformity (i.e. (3'). We may use then Eq. [3] with: 



From general considerations about the mechanisms generating the alignments one should 
expect this expression to be valid not only for the voids used to calibrate this expression 
but also for different voids, and also for halo alignments with respect to other large scale 
structures. 

Using Eq. [3] we may obtain any statistic concerning the orientation of the halos with 
respect to the voids. For example, the alignment of the minor and middle axes with the 
radial direction (that has been already measured in simulations) can be expressed in terms 




(3) 



0' 



1.22 (p- 1) 



(4) 
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of expression [3j The probability distribution, P min (cos7), for the cosine of the angle 7 
between the minor axis and the radial direction of the void is related to P or by the following 
expression: 

Pmin{cosi) = / / P or {cos 6, cos ip) 5 (sin 6 sin if) — cos ^d cos 8 d cos ip (5) 
Jo J 

where the Dirac's 5 distribution implement the constraint (based on spherical geometry) 
that only those values of 9 and ip corresponding to an orientation of the halo in which the 
cosine of the angle between the minor axis and the radial direction take the value cos 7 
contribute to the integral. P^. is simply the distribution implied by P or for the variables 9 
and ip: 

Carrying out the integral over ip we find: 



r^/ 2 , J sin 2 7 — cos 2 9 d9 

P min (cos7) = C0S7 / P or (cos9,^ — ) , — (7) 

A/2-7 sm9 A / s in 2 7-cos 2 ^ 



In a similar fashion, simply changing sinip by cosi(j in the argument of the Dirac 5 
distribution, we find for the probability distribution, P mi d(cos rj), for the cosine of the angle 
77 between the middle axis and the radial direction: 

rn/2 

Pmid(cosri) = / P or (cos9, cos 77/ sin 9)d9 (8) 

J-n/2-r) 

Using here expression [3] with =-0.134 and p=l.ll, which is the value corresponding to 
the halos studied by Cuesta et al. (2008) we recover (see Fig. |2]) the distributions P m id and 
Pmin determined directly in that work. 



3 ALIGNMENT OF THE ANGULAR MOMENTUM OF THE HALO 
WITH RESPECT TO THE HALO AXES 

The alignment of the angular momentum of the halo with respect to their axes have been 
studied in detail in simulations (Bailin & Steinmetz 2005; Gottlober & Turchaninov 2006; 
Patiri et al. 2006). We shall now consider a simple analytic model able to fit these alignments. 
To this end we take into account a set of halos with principal moments of inertia Ii, I2 and I3. 
Assuming that halos can interchange rotational kinetic energy and that they have reached 
statistical equilibrium, the probability distribution for the three components of the angular 
momentum J 1; J 2 and J3 along the principal axes is proportional to the Boltzmann factor: 
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Figure 2. Upper panel. The solid line is the probability distribution, P m i n (cos7), for the cosine of the angle 7 between the 
minor axis and the radial direction of the void as given by expression [5] with (3 =-0.134 and p=l.ll (see text for details). The 
dashed line is the analytic fit to the results from the simulations by Cuesta et al. (2008) given by Eq. 1 using p=0.884. Lower 
panel. The solid line is probability distribution, P m i,j(cos7), for the cosine of the angle 7 between the middle axis and the 
radial direction of the void as given by expression \E\ of /3 =-0.134 and p=l.ll (see text for details). The dashed line is the 
analytic fit to the results from the simulations by Cuesta et al. (2008) given by Eq. 1 using p=1.023. 



P(J)ocexp|-^^/fcT| (9) 

Let's write the components of J in terms of its modulus J, the angle 9j between J and the 
halo major axis, and the angle 0j between the projection of J onto the plane perpendicular 
to the halo major axis and the halo minor axis. Integrating over J we find for the probability 
distribution P j (cos 9 j,4>j)'- 



I, 1 + (/3- 1 - I, 1 ) cos 2 <f>j+(l- I, 1 - (/3- 1 - J^ 1 ) COS 2 <f>j) COS 2 9 
with J 2 = £ and J 3 = g and 9j e [0, vr/2] and 0j e [0, tt/2]. 



Pj (cos 9j,(pj) 

' (10) 
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If the distribution of rotational kinetic energy amongst halos were in statistical equi- 
librium, setting in this expression the mean values of 7 2 and I3 (72=1-38; 73=1.58, derived 
from the results of Jing & Suto 2000) one would obtain a good approximation to the actual 
distribution, since the dispersions of 7 2 and 7 3 are not too large. However, this is not what 
is observed in the simulations, because the actual distribution of the angular momentum 
is far from statistical equilibrium. Nonetheless, Eq. [10] can be used as a good model for 
the angular momentum distribution found in the simulations using the following values: 
7 2 =2.74 and 7 3 =7.41. On doing this, we find: < cos6»ji >=0.315; < cos6» J2 >=0.471 and 
< cos6*j3 >=0.671 where 6j\, 6*j 2 and 9j3 stands, respectively, for the angle between the 
angular momentum vector and the major, middle and minor halo axes. These numbers can 
be compared with the values found in the simulations for randomly chosen halos: 0.32, 0.48 
and 0.69 (Cuesta et al. 2008). In those simulations it has been also found that above 75% 
of halos have 9j\ larger than 60 degrees. Eq. [TO] predicts 73.5%. 



3.1 Alignment of the angular momentum of the halo with respect to the halo 
axes for a given orientation of the halo with respect to the void 

The above expression gives the probability distribution for the orientation of J with respect 
to the halo axes regardless of the orientation of the halo itself with respect to the large scale 
structure. However, from theoretical considerations we know that the distribution of 9j and 
cf)j depends on the orientation of the halo. We shall show later (in the next subsection) 
that for halos in the shell of the voids considered by Cuesta et al. (2008), the probability 
distribution Pj(9j,(pj\9,il ! ) for the angles 8j and <pj in those halos whose major axis form an 
angle 9 with the radial direction of the void may be approximated by: 



P J (cos9 J , ( Pj\9,ij) = -— 

71/ 2 

(7T x F2 x F3)- 1 / 2 

[0.365F2 + (0.135F3 - 0.365F2) cos 2 0j + (71 - 0.365F2 - (0.135F3 - 0.365F2) cos 2 fa) cos 2 9j} 3 / 2 

with Fl=l+a(cos 2 9-1 /3), F2=l+a(sin 2 #cos 2 /3), F3=l+a(sin 2 flsin 2 tp-1 /3) and a=- 
0.68 (9j G [0,7r/2] and <ftj G [0,7r/2]). This expression implies that for #=0 and any value 
of i[) (which corresponds to have the major axis oriented in the void radial direction with 
the other two axes lying in the void plane): < cos^ji >=0.251; < cos#j 2 >=0.489 and 
< cos 9j3 >=0.693. For 9=it/2 and ip—ir/2 (which corresponds to have the major and the 
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middle axes in the plane of the void and the minor axis oriented in the void radial direction) 
we find < cos #ji >=0.354; < cos 9 J2 >=0.526 and < cos#j 3 >=0.598. 

From the above numbers we see that the mean values of < cos 6^! >, being 9j\ the angle 
between J and the halo major axis, depends very substantially on the orientation of the halo 
with respect to the void. To understand the origin of this dependence we only need to consider 
the following: the axis corresponding to the smallest proper value of the shear of the velocity 
field associated with the void expansion is parallel to the radial direction of the void. This 
produces an extra component of J which is perpendicular to that direction. Consequently, 
when the halo major axis is parallel to the radial direction (6=0), the extra component of J is 
also perpendicular to the halo major axis. So, the average value < cos 9j\ > for this particular 
configuration must be smaller than the average value over all halos: < cos#ji >=0.315. On 
the other hand, when 9=ir/2, the tendency of J to lay parallel to the surface of the void 
does not imply an obvious tendency for < cos#ji >. But, to compensate for the smaller 
value of < cos#ji > when #=0, the mean value of < cos#ji > when 9=7r/2 must be larger 
than for a randomly chosen halo. 

The value of a given in expression [TT] corresponds to the halos considered by Cuesta et 
al. (2008), but that expression can be made universally valid by simply expressing a as a 
function of p (the parameter measuring the strength of the alignment of the major axis with 
the radial direction): 

a = -8.64(p - 1) (12) 

This relation should be valid in a general situation (with some reservations as presented in 
the discussion section) where both the orientation of the halos and their angular momentum 
are affected by a large scale shear of the velocity field. 

It must be noted that, unlike Eq. [101 expression [TT1 has not been checked with numerical 
simulations. It has been derived from Eq. [10] assuming a simple model for the dependencies 
on 9 and ip and fitting the value of a. In this respect it may be considered a prediction for 
the dependence on 9 and ip of the probability distribution for 9j and 4>j. 

3.2 Alignment of the angular momentum of the halo with respect to the void 
radial direction 

Having determined the probability distribution for the orientation of the halo with respect 
to the void, P or (cos 9, ip), and the probability distribution for the orientation of the angular 
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momentum with respect to the halo for a given orientation of the halo with respect to the 
void, Pj (cos 9j,4>j | 9, if)), we can now express the probability distribution Pj(cos/i) for the 
angle between J and the radial direction of the void in the following manner: 

Pj(coszi) = J- / / r r P or (cos9, 0)Pj(cos 9j,0j | 0,V) 
64 J -i J-i Jo Jo 

<5(cos 9 cos 9 j + sin # sin 9j cos(ip — 0j) — cos fj,)d<pjdipd cos 06? cos 9j (13) 

The factor 1/64 comes from the fact that the distributions P or and Pj are normalised 
with all angles between and ir/2. This is because, as long as we consider only the orientation 
of halos or the orientation of J within the halo, given the symmetries, this is the meaningful 
range for the angles. However, when we consider both orientations simultaneously to obtain 
the position of J with respect to the void we must consider all angles between and tt to 
cover all possible configurations. The factor 2 comes from normalising Pj (cos /x) to the range 
cos \x G[0,1]. 

Integrating over 0j and taking advantage of the symmetries of the integral, we find: 

— I f n /"l rcos(8—fi) 

P J (cos yU ) = -/ / / 

Jo J-l Jcos(9+fi) 

P or (cos 9, ip)P j (cos 9 j, (j) j | 9 ) ip)dcos9jdcos9dip 

i \ ) 

J (sin 9 sin 9j) 2 — (cos /i — cos 9 cos 9j) 2 

with <pj = arccos((cos/i — cos 9 cos 9j)j (sin 9 sin 9 j)) + ■0. 

If we assumed that the orientation of the angular momentum with respect to the halo 
is not affected by the orientation of the halo with respect to the void, the conditional 
distribution, Pj(cos#j,0j | 9,ip), would be equal to the unconditional one Pj (cos 6*j, 0j) 
given by Eq. [10] with I 2 =2.7A and 13=7.41. Then using Eq. [3] with /5 =-0.134 for P or , which 
corresponds to the halos around the large voids studied by Cuesta et al. (2008), and inserting 
both distributions into expression [13] we would not recover the distribution for cos /j found 
in the simulations, which is well parametrised by: 

Pj(cOS ^ )= (l + (^-l)cosV) 3 / 2 (15) 
with p=1.029. Instead, we obtain an alignment with the opposite sign (see Fig. [3]). 

To reproduce the observed distribution we must allow for a dependence of Pj on 9 and 

ip. A simple model for this dependence has been given in the previous subsection (see Eq. 

[TT]) . Inserting this into expression [T4l together with expression [3] with =-0.134 for P or we 
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Figure 3. The thick solid line is the probability distribution, Pj(cos/i), for the cosine of the angle fi between the angular 
momentum vector of the halo J and the radial direction of the void as given by Eg. 1141 with P or given by Eq.[3]and Pj given 
by Eq. 1111 The dot line represents the unconditional model which is estimated as above but using all F's equal to one in Eo llll 
The dashed line is the analytic fit to the results from the simulations by Cuesta et al. (2008) given by Eq. 1 using p=1.029. 

may fit a so as to obtain the best fit (with a=-0.68) to the numerical results which are well 
described by Eq. [151 This is how we have obtained the value of a given in [TIJ It must be 
noted that this expression represent a model in which, for any value of 9, ip the alignment of 
J with the axes follows expression [TU1 (corresponding to randomly chosen halos), but where 
the parameters show a small dependence on 9 and ip (the simplest with the right symmetry). 

The fact that expression UM with expression [51 and [Til for P or , Pj, respectively, agrees 
with [15] to a high degree of accuracy for any value of cos fi confirms the validity of this model 
for the conditional Pj (see Fig. [3]). 

If the two distributions entering expression [141 were independent of ip and i.e. if they 
take the following form: 



Pi(0,V) = --Pi(cos0) 

7T 

P 2 (9,?p) = -P 2 (cos0, 7 | 9) 



(16) 
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both normalised with all angles between and 7r, then expression [TH after integrating 
over ip, reduces to: 

1 f 1 /■«»(*-**) P 1 (cos9)P 2 (cos9 J I 9)d cos 9 jd cos 9 , . 

P(cos/i) = — / = (17) 

2n J -i Jcos(6»+ M ) y^sin 2 6> sin 2 9j - (cos fi-cos9 cos 9j) 2 

This expression gives the probability distribution for the angle [i between a fixed vector 
H and a vector a determined as follows: choose a vector b which forms and angle 9 with H 
with a probability distribution given by Pi and a randomly chosen (uniformly) azimuthal 
angle. Then, choose the vector a forming an angle 9j with b with a probability distribution 
given by P2 and a randomly (uniformly) distributed azimuthal angle. 

In the case where P2 does not depend on 9 this expression is symmetric in Pi and P2. 
This is not evident at first sight, given the different integration limits for 9 and 9j. P(cos/x) 
can be seen as the two-dimensional convolution on the sphere of Pi and P2. We will use 
this expression in the next section to study the effect of the dilution of the signal of the 
alignments due to: a) the miss-aligning processes after the moment of turn-around and b) 
the observational errors in the estimation of 9. 



4 DEGRADATION OF THE ALIGNMENTS 

4.1 Degradation of the alignments due to the miss-aligning process after the 
moment of turn-around 

In the standard (linear) tidal torque theory (TTT), the torque acting on a protohalo keeps 
its direction invariant through the time. It has been found in numerical simulations that the 
TTT predicts correctly the angular momentum of the protohalo up to the moment of the 
turn-around (Porciani et al. 2002). Hereafter, the tidal torque starts to rotate so by the time 
the halo virialises, its angular momentum will form an angle ( with respect to the initial 
torque. The initial torque direction is close to the direction of the angular momentum at 
the moment of the turn-around. This evolutionary process has been studied in detail by 
Porciani et al. (2002). We have found that the distribution of the angle ( that they found 
can be approximated by: 

P(C) ^ 16.83[e- sinC /°- 21 (l - 0(C - tt/2)) + 0.00856(C - tt/2)] sinC (18) 
where is the Heaviside function. 

Assume now that we have derived the theoretical predictions for the alignment of J with 
the large scale structure under the assumption of a non-rotating torque, whose probability 
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Figure 4. The thick solid line is the probability distribution, Pj(cos/i), for the cosine of the angle fi between the angular 
momentum vector of the halo J and the radial direction of the void accounting for the tidal torque rotation during and after the 
formation of the halo. To estimate the effect of the rotation of the tidal torque, we have assumed that, under the assumption 
of a non-rotating torque, Pj(cos/i) is given by Eg. 1151 with p=1.06 and we have convolved this expression with Eg. 1181 using 
Eq. 1171 The dashed line is the analytic fit to the results from the simulations by Cuesta et al. (2008) given by Eq. 1 using 
p=1.03. We find that the strength of the alignment, which is measured by p-1, reduces to a half with respect to the theoretical 
prediction under the assumption of a non-rotating torque. 



distribution we represent by P 3 (that turns out to be given by Eq. [15] with a certain value of 
p). To obtain the probability distribution, Pj, corresponding to the alignment of the actual 
J we must take into account the fact that the angular momentum of the halo is rotated by 
an angle ( with respect to the theoretical prediction. Consequently, Pj may be expressed as 
the convolution on the sphere of Pj and 2P(£), which is given by expression [T71 with these 
two distributions in the place of Pi and P2. The factor 2, in 2P(Q, is due to the fact that 
P(C), which is not symmetric with respect to tt/2, is not normalised within the interval 
(0,7r/2), unlike the distribution appearing in Eq. [T71 

If we use for P 3 a value of p=1.06, the resulting Pj (see Fig. HI) is also well approximated 
by expression [T5l with a value of p of 1.03. So, we see that the strength of the alignment, 
which is measured by p-1, reduces to a half with respect to the theoretical predictions with 
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a non rotating torque. This ratio between the value of p-1 for Pj and Pj is constant to a 
very good approximation for all interesting p values (p-l<Cl). 

Note that in obtaining Pj from Pj and P(() we have assumed that the rotation of J 
with respect to its primordial direction is independent of the orientation with respect to the 
large scale structure. This is a good approximation, because the large scale makes up only 
a small contribution to the total tidal field, and this contribution rotates less than that due 
to smaller scales. 

4.2 Degradation of the alignment due to observational errors at estimating 8 

When measuring alignments between the large scale structure and the rotational axis of 
spiral galaxies, the error in the estimation of the alignment angle due to uncertainties in the 
position of the galaxies and in the orientation of its axis leads also to a degradation of the 
intrinsic alignment. 

To quantify the relevance of this potential source of degradation of the alignment, we 
assume that the probability distribution, P er (Q, for the angle ( between the observed and 
the actual direction of J takes the form: 

Jq e^ 7 sin fidfi 

Assuming again that the distribution of the angle ( is independent of the large scale 
structure, we may insert in expression [T71 2P er (() and the probability distribution for the 
intrinsic alignment of the axes, P; n (cos/z), to obtain the distribution for the observed fi, 
P ;,(cos/i). Using for P in expression [T5l with p=1.06, we have obtained P Q b for four values 
of a (see Fig. \5§. Expected typical uncertainties in the observations, a ~20 degrees, do not 
produce a strong decrease in the signal. 

5 DISCUSSION AND CONCLUSIONS 

We have presented a set of general tools for the description of the alignments of halos and 
galaxies with respect to the large scale structure where they are embedded. Using these tools 
we have done the following: 

• We have provided an accurate model for the conditional probability distribution func- 
tion for the orientation of the angular momentum with respect to the halo axis for halos 
with a given orientation with respect to the large scale structure. We have fitted the single 
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Figure 5. The different lines shows the probability distribution, P b s (fJ.), for the angle fi between the observed angular 
momentum vector of the spiral galaxies J and the radial direction of the void assuming four different levels of uncertainty on 
the angular position of J (c=10; cr=15; <r=25 and <r=35 degrees). The distribution for the intrinsic alignment of the angular 
momentum is given by Eg. 1151 with p=1.06, whereas the observed distribution is characterised by the convolution, as provided 
by Eq. 1171 between the intrinsic signal and the error function given by Eq. 1191 



parameter of our model so as to obtain, using the appropriate mathematical identity, the 
alignment between the angular momentum of the halos in the shells of cosmic voids that 
has been found in simulations. The conditional PDF obtained in this manner shows a strong 
dependence on the orientation of the halo with respect to the void. This dependence is pro- 
duced by the tidal field on the shell of the voids which causes the minor axis to be more 
aligned with the radial direction than in the isotropic case, while the opposite is true for 
angular momentum. If that dependence did not exist, that is, if the PDF for the orientation 
of the angular momentum with respect to halo axis for halos in the shells of voids were the 
same as for randomly chosen halos (i.e. following an unconditional PDF) then the angular 
momentum would show the same kind of alignment with the radial direction than the mi- 
nor axis and opposite to that found in simulations. In the past, and under the assumption 
that the angular momentum PDF is independent of the environment, there has been some 
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confusion when interpreting the alignments of halo axis and angular momentum with the 
void radial direction. 

• We have shown that the strength of the predicted alignment of the angular momentum 
with the void radial direction, as given by the TTT, it is twice as large as the one actually 
measured in simulations (as parametrised in this work by p-1). To model this degradation 
of the signal we have convolved the angular momentum PDF characterised by a certain 
parameter p (Eq. [1]) with an analytical fit (based on previous results) that describes the PDF 
for the angle between the initial torque and the present angular momentum (Porciani et al. 
2002). As long as p-l<Cl the resulting distribution is also well described by the expression 
given in Eq. [1] with a value of p' such that l-p'=(l-p)/2. Understanding this degradation 
is a necessary step towards a fully theoretical explanation of the alignment of the angular 
momentum with the void radial direction found in simulations. In a future work (Betancort- 
Rijo & Trujillo 2010) we shall show that the TTT predicts, for the voids under consideration, 
a value of the strength of p~1.06 for the initial alignment, so a value of p~1.03 should 
correspond to the alignment of the present angular momentum, as found in simulations. 

• We have also convolved the initial angular momentum PDF with a Gaussian PDF to 
model the observational errors on measuring the angle between the angular momentum and 
the void radial direction. For real galaxies this error arises both from the uncertainty on 
estimating the direction of the galaxy rotation axis as well as from the error on determining 
the void radial direction due to the redshift distortion. We have shown that up to a r.m.s. 
of roughly 15% the intrinsic alignment is basically unaffected. 

In the present work we have dealt explicitly with large cosmic voids, but all expressions 
may be used whenever the orientation of the halos or galaxies with respect to the large 
scale structure (or more precisely, the orientation with respect to the proper axes of the 
deformation tensor or the tidal tensor on large scale) is considered. In order to use the 
expressions proposed in this paper in more general cases, two considerations are in order. 
First, the PDF for the orientation of halos (Eq. [3]) may show in general a dependence on 
<f). Therefore, in any expression involving integrals of P or , the dependence on <ft must be 
considered and that expression must also be integrated over and divided by 2n. Second, 
the ratios f3'/(p-l) (see expression HI) and a/ (p-1) (see expression [T2l) have been determined 
using halos within certain shells of cosmic voids. However, from general considerations about 
the mechanisms generating the alignments it follows that these ratios should be essentially 
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independent of the voids used (as long as p-l<0). In other cases (i.e. walls or filaments) the 
same considerations lead to the constancy of the first ratio. The effect described by the second 
ratio (i.e. the dependence on environment of the orientation of the angular momentum with 
respect to the halo axes) does not depend just on the present deformation tensor, but also 
on its deformation history (we will show that on Betancort-Rijo & Trujillo 2010). Therefore, 
its value could depend on the kind of structure considered. This dependence, though, should 
be expected to be small, but as much as it can be measured it might provide an interesting 
test for the standard (and non-standards) model of the large-scale structure formation. 
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